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' Abstract: We consider a class of mechanical particle systems with determinis- 

tic particle-disk interactions coupled to Gibbs heat reservoirs at possibly differ- 
<~i ' ent temperatures. We show that there exists a unique (non-equilibrium) steady 

OhI state. This steady state is mixing, but not exponentially mixing, and all initial 

distributions converge to it. In addition, for a class of initial distributions, the 
rates of converge to the steady state are sub-exponential. 



1. Introduction 



Rigorous derivations of heat conduction laws for mechanical particle models 
, coupled to heat reservoirs remain a mathematical challenge. A variety of models 

QQ ■ have been introduced in the past [TI^M fTTlfTSlI 1 6 lIlT] ; nearly all of the proposed 

I derivations of the Fourier Law are partial solutions based on unproven assump- 

■ tions [2l[3ll4l[5l [m[T8] . Developing proofs of these assumptions would require deep 

' understanding of the properties of systems in non-equilibrium, i.e., coupled to 

, several unequal heat reservoirs. The standard assumptions include the existence 

ff^ ' of the unique invariant measure (steady state) as well as certain bounds on the 

rates of convergence of initial distributions to the invariant measure. 

For systems in equilibrium, i.e., when the temperatures of all the reservoirs 
are the same, the steady states can often be written down explicitly. The question 
. of existence of non-equilibrium steady states has been open for practically all 

^ ' mechanical particle systems, by which we mean Hamiltonian-like systems, driven 

- - - by stochastic heat reservoirs. The main difficulty lies in dealing with the non- 

compactness of the phase space. For the systems under consideration, however, it 
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is relatively easy to envision scenarios under which particles slow down (freezing) 
or speed up (heating) due to stochasticity of the heat reservoirs. This may push 
initial distributions towards zero or infinite energy levels and ultimately violate 
existence of physically relevant steady states. 

An example of freezing has been observed (numerically) in one of the proposed 
models [12 . In that model, a particle acquires very low values of the speed once in 
a while under the evolution of the dynamics (due to stochasticity). Low values 
of the speed lead to long traveling times between collisions during which the 
particle has no influence on the evolution of the system. It is observed for the 
system in JT2J that more and more particles get stuck on the low energy states 
resulting in fewer and fewer collisions per unit time. To rule out such unfortunate 
scenarios one must be able to control the probabilities of particles acquiring low 
speeds and the rates at which the speeds recover to normal ranges. 

The dynamics of the mechanical particle systems driven by heat reservoirs 
may be viewed as a continuous-time Markov process. Harris' ergodic theorem 
and its generalizations are common tools for obtaining existence and uniqueness 
of the steady states as well as (exponential) convergence of initial distributions 
to the steady state. In discrete time, the theorem requires two things: to produce 
a non-negative potential V (Lyapunov function) on the phase space which, on 
average, decreases as a power law under the push forwards of the dynamics, and, 
given such a y, to show minorization or Doeblin's condition on certain level 
set of V. The first condition guarantees that the dynamics enters the 'center' 
of the phase space, a certain level sets of V, with good control on the rates; 
and, once at the 'center', coupling is guaranteed by the minorization condition. 
Harris' ergodic theorem was applied in (14, for a discretization of the original 
continuous-time Markov process. In j23j the results were extended to continuous 
time. The existence of a unique steady state was obtained by constructing a 
suspension flow over the discrete dynamics; convergence of initial distributions 
to the steady state followed from a result in general state space Markov process 
theory |21) once irreducibility of time-1 discrete process was shown. In addition, 
P5] demonstrated that this convergence is sub-exponential for a class of initial 
distributions. The slow rates of convergence are due to the abundance of slow 
particles in the system which, in turn, do not influence the system for extended 
periods of time. This slows the rates of mixing. 

The analysis in [Hll^ relies heavily on the fact that there exists a mean- 
ingful discretization of the system that mixes exponentially fast. Because parti- 
cles do not interact, the study of the dynamics of one particle on the collision 
manifold reveals important dynamical properties that yield implications for the 
continuous-time system. For an interacting particle system, the rates of mix- 
ing for the continuous-time system and its relevant discretizations happen to be 
comparable due to the slow particle effect. 

Sub-exponential mixing seems to be prevalent for canonical interacting par- 
ticle systems driven by Gibbs heat reservoirs. If a discrete system does not mix 
exponentially, finding the potential V that would still guarantee existence of 
invariant measures is a very diligent task requiring extremely good understand- 
ing of the dynamics of the system. Thus, one needs different methods to tackle 
the question of existence of non-equilibrium steady states. The paper of Meyn 
and Tweedie [19) provides a general framework of showing existence of invariant 
probability measures for general state space Markov processes. 
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In this paper we consider a class of mechanical systems in which particles 
interact with an 'energy tank' represented by a rotating disk anchored at the 
center. Particles move freely between collisions with the tank or the boundary. 
When a particle collides with the disk, an energy exchange occurs, in which the 
particle exchanges the tangential component of its velocity with the angular ve- 
locity of the disk and the normal component of the particle's velocity changes 
sign. A system in this class is coupled to heat reservoirs set at possibly differ- 
ent temperatures that absorb particles when they collide with the boundaries 
of the reservoirs and emit new particles according to the Gibbs distribution 
corresponding to the temperatures of the reservoirs. 

The main geometric assumption that makes the analysis feasible is that a 
particle can hit the disk at most once before returning to the heat reservoir. The 
existence is shown in section [3] through estimating hitting times of a carefully 
chosen compact set C in continuous time without an aid of a discretization or 
a potential. A regeneration times idea is employed in the argument. In order to 
apply a general state space Markov process theory developed in [TU], one also 
needs to show that the minorization condition holds on C, which we do in 13.41 
Convergence of initial distributions to the steady state follows by application 
of a theorem in |21| after a small modification of the minorization condition 
argument. In section 0] we show that mixing is not exponential and for a large 
class of initial distribution convergence of initial distributions to the steady state 
occurs at sub-exponential rates. The argument is similar to |23| . However, we 
cannot use the potential V for certain upper bound estimates and different 
methods are required. The key property that leads to sub-exponential mixing 
is that the (invariant) measure of the states for which at least one particle will 
not collide with a heat reservoir or a disk in time t is of the order of r~^. 
The dynamics thus resembles one of an expanding map with a neutral fixed 
point and sub-exponential convergence rates for a class of initial distributions 
can be obtained using arguments similar to |24j . The bounds on the measure of 
the particles that will not collide with a heat reservoir or a disk in time r are 
obtained using the minorization condition and other properties of the dynamics. 



2. Settings 

Let 7^ be a circular domain of radius R + d. A disk D is anchored at the center 
of dr. The disk D is allowed to rotate freely with angular velocity cj; denote 
by e S*^ a position of a marked point on dD. Let i„ and Ld be two vertical 
walls on the top and on the bottom of D, splitting F \ D into two halves, 



left and right. See Fig. 1(a) Our system consists of k particles m. F \ D with 
positions {xi,X2, • ■ ■ Xk) =■ x and velocities (vi, f2, ■ • • , Vk) =■ v. The particles 
are confined to either of the two halves and move freely between collisions with 
dF U dD U L„ U Ld- The collisions with L„ and Ld are specular, i.e., the angles 
of incidence are equal to the angles of reflection. When a particle collides with 
the boundary of the disk dD, an energy exchange occurs, in which the particle 
exchanges the tangential component of its velocity with the angular velocity of 
the disk and the normal component of the particle's velocity changes sign. More 
precisely, if u = ivt,v±) is the particle's velocity decomposition upon collision 
with dD and the disk rotates with angular velocity w, then the post-collision 
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(a) The configuration (b) Suspension flow coordinates 

Fig. 1. The system 



velocities are: 

v'^ — — w_L, v[ — Lu, and uj' = vt, 

where v' = is the particles velocity decomposition and uj' is the disk's 

angular velocity immediately after the collision. This interaction was introduced 
in \15Mi- 

The left and right parts of dF, dFi and 9I2, act as heat baths at possi- 
bly different temperatures Ti and T2 respectively. Particles get absorbed by 
the heat baths upon collision with dF, and, upon collision of a particle with 
the disk, a new particle is emitted immediately at the collision location with 
speed s S (0, 00) and angle (p S (— fif) distributed according to d{s,(p) — 

s^e~^j'')^ cos{ip)dsdip, where = 1 or 2 depending on whether the parti- 
cle is confined to the left or the right half of F. This distribution on the boundary 
dF corresponds to particle's velocity distributed as ^i^e~'^^('' dv in F. 

We would like to define the associated Markov process with the dynamical 
rules governed as above. A phase space for such a process should consist of 
quadruples (x, v, 6, w) with proper identifications at the collisions. In particular, 
when Xi € U i^, if has positive (negative) horizontal component, then the 
corresponding particle is confined to the right (left) half of the domain F\D. To 
simplify the notation, we would like to exclude states {(x, v,0,a;) : for some i, 
Xi € LuULd and Vi has zero horizontal component} from the phase space. In our 
future arguments we will frequently omit mentioning the situation of particles 
reflecting from U Ld'. due to symmetry and circular shape of F, the distance 
of flight is the same whether reflection occurs or not. 

In addition, we would like to exclude all the states with stopped particles 
{vi = for some i) and all the states that lead to such with positive probability. 
A particle stops if it collides tangentially with a stopped disk {lu — 0). Conse- 
quently, any state (x, v,0,w) with a particle heading for a tangential collision 
may lead to positive probability of stopping a particle as time evolves. Let 

f2 ={(x,v,6',w) :v,^0,l<i<k, 

if e L„ U Fd, then the horizontal component of Vi is nonzero, and 
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no particle is heading for a tangential collision with the disk} / ^ 

where ~ corresponds to a choice of outgoing velocities upon collision of a particle 
with dru dDU LuU L^. 

Let <?t be the associated Markov process on i7; denote the transition proba- 
bility kernel by V* . Note that SI is locally compact and separable, and <l>t has 
strongly Markov with right-continuous sample paths because we chose to keep 
track of the outgoing velocities are collisions. Those assumptions are necessary 
in order to apply general state Markov process theory in section [31 

Theorem 1. There exists a unique invariant measure for the Markov Process 
<Pt ■ It is mixing, hut not exponentially mixing. Moreover, all initial distributions 
converge to /i, but for a large class of initial distributions the convergence rate 
is at best polynomial. 

In section[3]we will show existence of an invariant probability measure, which 
is mixing; in section 2] that mixing and the convergence of initial distributions 
to /i for a class of initial distributions are not exponential. 

3. Existence and Mixing 

Definition 1. A non-empty compact set C G B{f2) is called mc -petite if there 
exist T,ri > such that 

V^ix, ■) > r]mc. 
Here mc is the uniform probability measure on C . 

The condition V"^{x, •) > rjmc is frequently called the minorization or Doe- 
blin's condition on C. 

For any 6 > and a set C G B{n) define tc to be the first hitting time on C 
and Tc{S) to be the first hitting time on C after waiting time 6. 

We will use the following result by Meyn and Tweedie on continuous-time 
general state Markov chains [H] : 

Theorem 2 ([19]). Assume there exists a petite set C such that 

^ Pz{tc < oo} = 1 for all z £ f2 and 
— for some S > 0, supEz[Tc((5)] < oo 
zee 

Then there exists an invariant probability measure for <l>t. 

We obtain bounds for the expected values of the hitting time of C using 
the regeneration times idea explained in 13.11 Roughly speaking, with the aid of 
stopping times, we split each random trajectory into similarly behaving pieces 
of bounded, on average, length. We define the set C and the stopping time r 
and show how regeneration times idea applies in our setting in 13.21 We estimate 
the expected times of the lengths of the pieces in l3.3l In l3.4l we show that the set 
C we defined is indeed petite. This ensures the existence of a unique invariant 
measure fi. 

We show that fj, is mixing and all initial distributions converge to it in 13.51 
using a theorem |211 Thm. 6.1] of Meyn and Tweedie that applies to continuous 
times Markov processes on general state spaces. To satisfy the conditions on 
the theorem one needs to demonstrate that some skeleton chain of the Markov 
process is irreducible. We show that the time-1 chain <Pi is irreducible using 
a modification of the argument for showing that C is petite. 
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3.1. Regeneration times idea. To perform the estimates required for Theorem 
m we apply the idea of regeneration times. We wih first describe an ideahzed 
situation, and then show that a very similar argument works for the real 
Assume for a moment that there exists a distribution v and a stopping time r 
such that {<l>t+Tit > 0} is independent from > 0} as well as t > 0} 

is stochastically equivalent to {^t,t > 0} in a sense that they have the same 
joint distributions provided that is distributed with v. Then the process can 
be split into independent regeneration epochs {t„}. If, in addition, ¥,^[t\ < oo 
and a is a geometric random variable with p = v{C), then 



Theorem [5] asks for the initial distribution of to be a point measure at 
z e J7; assume there exists a stopping time f such that $f is distributed with v. 
If 

sup2,g(7{Ez[7^(<^)]} < 00 and Ejy[T] < oo, then by a similar estimate we conclude 
that sup2;g(7{IEz[Tc((5)] < oo}. If, in addition, we show that Pz[f < cxi] = 1 for 
all z e i7, then Pz {tc < oo} = 1 for all z e ^2 follows too. 

For the actual process ^t, the first part of the argument applies if v is the 
invariant measure. Not only we do not know if it exists, showing that there exists 
a stopping time f with <Pf distributed as v is nontrivial. However, if there exists 
a stopping time r such that the system 'almost renews' at time t since 'enough' 
initial data is forgotten by time r due to the randomness of the heat baths, the 
argument may still carry through. 

A bit more precisely we would like to find a stopping time r such that 
{^t+T,i > 0} is 'independent enough' from {^t,t > 0} and 'Pr is distributed 
'similar' to v given that is distributed 'similar' to v. This will guarantee al- 
most geometric rates of hitting C . In addition, we want r to be 'small enough' 

so that SUP2gc;{II^z['''('^)]} < °° ^^"i Ey[T] < CX). 

3.2. Proof of existence. Let z be a state in J?. Some of the particles in z may be 
heading for a collision with the disk in a sense that each of these particles will 
collide with the disk before colliding with dF. Let to(z) be the time of the last of 
those collisions with the disk given z. Note that to(z) is finite and deterministic. 

Let T be the minimum time at which all particles and the disk randomize. 
More precisely, starting with z e J?, r = min{i > : such that both of the 
following events have occurred: 

— all particles in z have collided with dF at least once (ensures that all the 
initial particles' velocities are forgotten); 

— a particle originated from dF, hit the disk at some time t > io(z), and collided 
with dF again (ensures that the angular velocity uj is forgotten) }. 

A priori it is not completely clear whether r is almost surely finite; we will 
show it along with the expected value estimates. Though, at time r, the initial 
velocities of the particles and the angular velocity of the disk are forgotten, the 
positions x and 9 may still be strongly correlated since collision times may be. 



= P[(T = 1]E^[ti] + ¥[<7 = 2]E^[ti + ra] + • • • 
v{C)%Ar] + (1 - v{C))i^{C)2¥.,[t] + ■■■ = 5^ 




< oo. 
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Also note that (Pr belongs to the collision manifold J7o = {z G i7 : a;^ G dF for 
some 1 <i <k}. 

In the following it is convenient to introduce a change of variables in order 
to make passing from to fio easier. We would like to replace (x, v, 6, to) with 
(r, s, iy9, 5, 0, w), where the new coordinates are based on the information from 
the past or the future collision. See Fig 1(b) More precisely, let r = (ri, • • • , r^), 



s = (si,---,Sfc), if = {(fi,- ■ ■ ,(fk), and ^ = (Ci, be as follows: is the 

point of the past collision of the particle with dF if its previous collision 
was with dF and not with dD (here we do not count collisions with L„ U Ld)] 
otherwise is the point of the future collision of the i*'* particle with dF. Note 
that the geometry is chosen in such a way that a particle can experience at most 
one collision with dD between collisions with dF. In the first scenario, > is 
the distance of flight of the particle to dF in the direction of —Vi (with possible 
reflection off L^U Ld, and in the second scenario £,i < is the distance of flight 
of the particle to dF in the direction of Vi. Let (p ^ (— ^, ^) be the angle with 
respect to the normal to dF at the collision point, and Si is the speed of the 
particle. We will denote states in !7 by z. We will also denote angles of collision 
with the disk by ip' . Note that sm{ip') = -^^^ sm{ip). 
Then 

SIq — [{(r, s, (p,£^,uj,d) ^ n : — for some i, 1 < i < k}. 

Suppose we start with a point measure Sz, z G f2, and run <Pt until time t. 
Once time r is reached, the speeds s; and angles 0i of the particles are distributed 
independently with 

d{sl,ipj) = pi3.^^^{si,Lpi)dsidip, := {2/3^1^^^/ y/Tr)sfe~'^^w''' cos{ipi)dsidipi if C» > 

and is distributed with d\uj'^\ = pfi.{uj)duj — (2-y/^/-\/7r)e~^^'^ doj, where 
j = 1 or 2 depending on the side from which the disk experienced its last 
coUision before time t (in case of ^ < 0, i.e. after disk colhsion, d{sj , ipj) is the 
same if = /3j and is a certain mix distribution if Pj). The inverse 

temperature Pj is the only memory kept for the distribution of \uj\ at time r. 

Let /iz := VJSz. We would like to show that due to 'enough randomization' 
of speeds and angular velocities, the expected time E^^ [r] has a uniform upper 
bound for all initial z. The same holds for all subsequent r-stops. More generally. 

Lemma 1. There exists D > such that for any initial distribution X, if v = 
Vl\, then^^lr] < D. 

Note that Lemma [T] does not say anything about the initial waiting time for 
renewal r starting from z S J7 or an arbitrary initial distribution A, only about 
the subsequent ones. We will prove Lemma [1] in subsection 13.31 

In order to apply Theorem [21 we need to find a petite set C such that proba- 
bilities of hitting C at regeneration times t„ are uniformly bounded away from 
and 1. In addition, we want sup2.g(7{Ez[T((5)]} < oo. Roughly speaking, we 
expect C to be a collection of all states z in with uniform upper and lower 
bounds on particles' speeds and an upper bound on disk's angular velocity both 
before and after collision. Since (1>t has values on Hq, but, in reality, we are 
concerned with C defined on whole fi, it would be convenient if C was forward 
invariant during particle's flight from dF to dF. 
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If a particle collides with a nearly stopped disk, w 0, nearly tangentially, 
Lp' ~ ±-|, then the particle's velocity is significantly reduced, s' « 0, and it 
would take a very long time for such a particle to reach dF . So when defining 
C, we would like to restrict disk collisions from being too close to tangential. In 
fact, we would like | sin((/j)| to be bounded away from -^==5., where a — 
this will be used later in the proof. 

Let C to be the set of all z G i7 such that for some Smin < Smax and e > 0, 
|w| < Sinaxy^i^- and for all i, 1 < i < fc: 

~ either | sin((p)| > a (no collision with the disk) and Smin < Si < Smax, 

— or I sin(iy9)| < --^===vT^^, > (before colhsion with the disk), and Smin < 

- or |sin(v3)| < . , , , =, 6 < (after collision with the 

disk), and Smin-y^a^^ — ■^i — Smax-y^Q,^^^^ • 

Proposition 1. C is petite. 

We will prove Prop. [1] in subsection 13.41 

Easy computation ensures that C is forward invariant between collisions of 
particles with dF. Moreover, the only way to leave C is for some particle to 
collide with dF and to originate with new speed and angle not satisfying the 
speed and angle conditions above. 

Similarly, in order for <Pr G C, all the particles' speeds and angles must be 
drawn in the correct range. Since the drawings are independent, there are a lower 
bound 7min and an upper bound 7rnax on the probability for <?t- G C. 

To show that the conditions of Theorem [5] are satisfied, we need to estimate 
the expected value of the first hitting time r provided we originate in C and 
wait some small time S. 

Lemma 2. There exists 6 > and D' > such that 

sup{E,[r(5)]} <D'. 

Moreover, Pz{r((5) < 00} = 1. 

We will prove Lemma [5] in subsection 13.31 

Let (7(7 be the hitting time of the set C for the stopped Markov chain <l>r , i.e. 
<Tc = mm{n > : {<PrT G C}. Then for z G C, 



^z[Tcm < P(<Tc = 1)E,[t(^)] +P(ac = 2){E^[T{5)]+Er:sAr]) 
+P{ac = 3)(E,[r((5)] + Ep.^^ [r] + Epj.^^ [r]) + • • • 

< 7maxi?' + (1 - 7mi„)7max(i^' + D) + {I ~ 7min)'7max(i?' + 2D) + ■ ■ ■ 

= ImaxD' [1 + (1 - 7min) + (1 - 7min)^ • • •] + (1 - 7min)7max£' [1 + 2(1 - 7„iin) H ] 

7niax ^/ ^ (1 7min)7max ^ 

7min 7min 

This is the second condition of the Theorem [2] The first condition follows 
from Pz{t{5) < 00} = 1 in Lemma[2]and the estimates above. We conclude that 
there exists an invariant probability measure for <l>f O 
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3.3. Estimates for Ez(t). Before proceeding with estimating Ez(t), let us first 
estimate the expected value of a flight time from dF to OF by a particle emitted 
randomly. Let T^^'s^"* be the random stopping time of hitting dF by a particle 
emitted at time with velocity distribution d{s, ip) = pj^. (s, (p)dsdip that will hit 

the disk rotating at angular velocity w. Let I := + d)^ — i?^ be half of the 
maximal distance of flight. Then 

/>oo 



sin((^)>a 

f f°° I I 

J Jo s ^/UJ^ + s-^ cos^[ip') 



Note that 



sin(c^)<a 



< 



I 



V^^ + s^cos^^') scos(^') s^l-sm\^)/a^' 
Therefore 

J Jo S ^/-K 

sin((/5)>a 

/ [-+ ]-^s^e-^^-" cos{^)dipds < K, 

sinM<a^° " Vl-sin^(^)/a2 

for some K > independent on lo. 

This guarantees, in particular, that every randomly emitted particle has finite 
expected time of flight. 

Let T^'* be a random stopping time of hitting dF after one disk collision by a 
particle emitted from dF with (s, (p) ^ pp. (s, (p)dsdip. The probability of hitting 
the disk at each round is a, so 

E[T^^] <aK+{l- a)a2K + (1 - afaSK + ■■■ 

= aK[l + 2(1 - a) + 3{1 - a^) + ■ ■ ■] = ^ = -. 

Now we are ready for preliminary estimates for Ez(r) for any z & f^. 

First note that r < maxi<i<fe{r*} + T^'s'^* + T'^'*, where T' the random time 
it takes for particle z in z to reach dF. Indeed, if no random particle collided 
with a disk before the last one hit dF, then either the last particle will hit the 
disk in time T^^^ or some other one will before; it will take less than T^'^ht ^img 
for it to exit. Therefore 

EzM < max {E^T']} + E[T^'^^^] + E[T^'^] 

l<i<k 
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Let us estimate Ez[T^]. The initial position and velocity of the particle uniquely 
determine whether the particle will collide with the disk or not. In case of no 
collision, the time of flight is < If the particle is headed for a collision with 
the disk, there are three possibilities for the value of the angular velocity of the 
disk upon collision: 

1. original w; 

2. SiSm{(p[), if the z*'' particle that collided with the disk immediately before 
the first particle; and 

3. a random angular velocity uj acquired by the disk due to collision with a 
particle emitted randomly from dr. 

In the 3'"'^ situation, \uj\ is drawn from the distribution 2\J^e~^^^ duj. The 
expected value for the exit time after collision with the disk is bounded as follows: 



/ / / B.j a r,2 I I 



W V TT VAninTT 

For the convenience of notation below, set 1^9' = if sin((/3) > j^^- Then 
E,[Ti] < l+max{-, ^ , ^ ^ ^ } 

^ 21 I 

Similar estimate holds for Ez[r'], 2 < i < k. 
Therefore 

r 1 r 2Z , I (l + a)K 

r < max { —} + _ + ^ '— 1 

l<-*<fe's,COs(.^^)' VAninTT 

Proof of Lemma [TJ 

The bound Ez[t] in the equation ([T|) depends only on Si and cos{ip'^), i — 
1,2. In the distribution = VI X, each Si and ipi are distributed according to 
Pp {i){si,'^i)dsidipi (before disk collision). 
'Then 



E, T < / max{ } + ^^ + \ >—]di, 

J 'l<*<fc's,C0s((/3^)-' V/^minTT a 

n 



< 



21 , , I {l + a)K 

-T:Pf},{s,^)drdLp ' 



SC0S((/?') ' ' V/^minTT 



arcsiii Q />oo ^7 OR'^^'^ 

— s e '^^ cos{Lp)dipds 



-^0 sy^l ~sin2((^)/a2 ^ 
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j^s'^e P''' cos{ip)dipds + + — < D 



J arcsin a J 

for some D > Q. O 



1^ ' ' VK 



Proof of Lemma [2j 

If z e C then 

^ r 1 r ^ . 2Z (l + a)K 

T < max { — -} + -j== + ^ '— 2 

21 I il + a)K 
< —p V —;= + — 

V/^minTT a 

If we start with z €E C and wait for some time 5, then some of the particles may 
experience coUisions with dF and redistribute their speeds and angles according 

to pp-(i){s,(p)drdip. Then 

2/ 

Ez[r(<5)] < ^ + D=:D', 

where D is the constant from Lemma [1] 

Now suppose we start with any z g ^2. The stopping time r < maxi<i<fc{T*}-f 
yfiight _j_2.hit_ rpj^g ^ijjjg j^z < Ji^^^^ + ^ is finite almost surely; so are E[T^'s^^] 

and E[T^'shtj gi^ce their expectations are bounded. For almost every y in the 
support of V^lz,-), T is also finite almost surely. Therefore, for any z € 12, 
Pz{t((5) < cx)} = 1. □ 



3.4- C is petite: proof of Prop. [J^ We would like to show that there exits T, 77 > 
such that for any z € C, 7^^(z, •) > rpnc, where mc is the uniform probability 
measure on C. This statement implies, in particular, that for any z, z' G C, 
there is a sample path a{z,z',T) which takes precisely time T to complete. We 
will start by showing this implication, with additional restrictions for the path 
to be 'regular' in a sense that it stays away from tangential collisions (precise 
definition to follow) and for a{z, z' ,T) C C. Those restrictions will be important 
later when we extend our argument to showing that C is petite. Our proof follows 
general outline of the proof of the minorization condition, Prop. 2 in |14| . 

The state z e C represents positions and velocities of k particles as well 
as the marked position and the angular velocity of the disk. Particles and the 
disk interact as the dynamics evolves. In order to make the analysis simpler, we 
would like to choose a path that 'decouples' particles and the disk. We achieve 
it by imposing a rule that as soon as a particle reaches OF, it does not collide 
with the disk anymore; the exceptions are that some particle needs to reset the 
disk's angular velocity and the final trip to z'. Note that, with this assumption, if 
z E C, the times ti for i*^ particle to reach dF for the first time are deterministic 
and uniformly bounded by some t > 0. We treat the final state z' similarly by 
running the dynamics backwards in time. Let t[ be the times for the particles in 
z' to reach OF when running the system backwards in time; note that tg < t. 

In addition, we impose that, at each collision with dF except setting the disk 
and the last one for each particle, the emission angles are uniformly bounded 
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away from and 1. Namely, we require that the emission angles satisfy a+e < 

sin{ip) < 1 — e for a — -j^^ and some e > 0. 

The 'decoupling' of the particles reduces the problem to a sub-problem con- 
cerning only one particle, which can be stated as follows: 

Lemma 3. For large enough T and for any r,r' € dF^ (or dF^), there exists a 
particle path Ui from r to r' , with outgoing angles and speeds satisfying a + e < 
|sin((^)| < 1 — e and Smin < s < Smax, such that Ui takes precisely time T to 
complete. Moreover, the paths can he chosen in such a way that the number of 
collisions is hounded by some monotone function of T . 

Proof of Lemma [3j 

Let TO be the minimal number of collisions required to travel between two 
diametrically opposite points satisfying the angle assumption a -I- e < | sm{Lp) \ < 
1 — e (usually m — 2, but if d is very small compared to R, m may be larger). 
Then to collisions is enough to travel between any r and r'. Let tmin and fmax 
be the minimal and maximal times of travel along such an path with < to 
collisions satisfying the speed and angle parameters of Lemma [31 Note that a 
simple application of the Lagrange multipliers method guarantees that, for all 
intermediate values t G [tmin, ^max], there exists a path taking precisely time t. 

Let iniin and imax bc the minimal and maximal times of travel along a back 
and forth path originating and ending at r', satisfying a -I- e < | sin(iy9)| < 1 — e 
(^min and imax do uot depend on the choice of r'). By appending this path to 
end of the path between r and r', we can produce paths from r to r' taking any 

t G [^min + fc^minj^max + fc^max] with TO -|- 2fc Segments. Let fc = miujfc : iniin < 

(imax - tmin) + ^(tmax " Wm)}- Then for any T > tmin + Hinin, there exists a 
path between r and r' taking time T to complete with the maximal number of 
segments to -I- 2[ -^r*""° J , as desired. □ 

Given Lemma [31 in order to reach z', we just need to send a particle to reset 
the disk and to reconnect the paths of particles matching them in an arbitrary 
manner in dF^ and in dFn. For particles that do not reset the disk, T = T—ti—t[. 

The particle that resets the disk has to change the angular velocity from a; 
to oj' , where w and w' are the angular velocities of the disk after all the particles 
in z or z' have collided with dF in forward and backward times respectively. To 
change w to uj' one simply needs to emit a particle with (r, ^, s, yj) G C such that 
w' = ssm{ip') — ssin((^)/Q;. However, for the purposes of obtaining the lower 
bounds in our subsequent argument, we would do it in two steps, by first setting 
to some intermediate value Cj and only then to uj' . In most situation we would 
choose w = 0. 

Let t^'^^ denote the time to reset the disk. Then < t^°^ < i!,'?' , for some 
*min 5 *max > • ^^r the particle that resets the disk we choose f ^T-t^-t'^- 1"*"* . 

If we take F > 2i + i^^^ + t,nin + fctmin, then we are guaranteed to have a 
path from z to z' with < TV := m -|- 2[ ''"^*'"'" J -|- 2 collisions with dF. Call such 
a path (^(z, z', T). 

Thus, we obtain: 

Lemma 4. Given e > 0, there exists F > and N > such that for any 
z, z' £ C, there exists a sample path cr(z,z',r) making less than N collisions 
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such that between the first and the last collisions of each particle with the dF , 
all Sniin < Si < Smax md a + € <\ sm{if)\ < 1 ~ £• 

Density bounds along f7(z,z',T). 

The next step is to sliow tliat if we start with, a point measure at z, S^,, 
it 'acquires density' as it evolves with the dynamics in a sense that VlS^ has 
a nontrivial absolutely continuous component for large enough t. Moreover, the 
density of this absolutely continuous component is uniformly bounded away from 
zero in a neighborhood of each path (t(z, z', T) and, in particular, at the endpoint 
z' G C. 

Density in a neighborhood of a point z G is the product of the densities 
in neighborhoods of coordinates of each particle and the disk. For the majority 
of the path (t(z, z',T), particles and the disk do not interact, so we can deal 
with their densities separately. In Lemma [5] we show that each particle acquires 
density with a uniform lower bound in a fixed size neighborhood of the second 
collision and LemmalHlkeeps track of lower bounds of the densities at subsequent 
collisions along the path as we push the measure forward. Lemmas [7] and |H] deal 
with the particle that collides with the disk, making the disk acquire density 
and the particle to re-acquire density at the next collision with dF after loosing 
some to the disk. Since the number of collisions along a{2,,-z! ,T) is bounded, if 
we combine Lemmas [S]|S] together, by the last collision, each particle as well as 
the disk have density in a fixed-size neighborhood of the collision point. After 
the last collision the dynamics is deterministic and the value of the density is 
preserved under push forwards. This allows us to conclude that VjSz has a 
uniform lower bound on the density at z'. 

To formalize the argument, we need to define more precisely what we mean by 
a neighborhood of a collision point: in the coordinate system defined in section 
[5] the dynamics has discontinuities at collision points. Let us first replace co- 
ordinate by sin((/j): this coordinate change makes sin((/j) re-distribute uniformly 
at collisions; in addition, for (r, sin(i^), ^)-coordinates, the Jacobian for the stan- 
dard billiard flow is equal to 1 (see for details) . Then, at each collision point 
of ct(z, z',T), we extend the coordinates (r, ^, s, sin((^)) forward or backward in 
time to accommodate neighborhoods of fixed size at each collision point for 
some small enough > 0. Here we assume that ^ is taking negative values before 
collisions and positive values after collisions. This extension is possible due to 
the bounds on the speeds and angles introduced in Lemma |4l 

First, we show that by pushing (5z forward, one can acquire density with a 
uniform lower bound. By the design of our path cr(z, z', T), as soon as a particle 
reaches dF ^ it is independent from other particles. Therefore, we only need to 
show that the each particle acquires density with a uniform lower bound in a 
^-neighborhood of some point along a particle sub-path in (7(z,z',T). 

Given z = (r, 0, s, sin(93)) G i7o and C > 0, let 

H^{r, 0, s, sin((^)) := (r-C, r+Q) x +C) x (s-C, ^ + x (sin((p)-C, sin(^)-l-C) 
and let ^^H^{rfl.s,s\n{^p)) be the uniform measure with density 1 on H. 

Lemma 5. There exist C^Vo such that F[i\Pl^5z\ > 'ilol^Hf{r.o,s.sin{ip)): 1 < * < 
where Ti be the time of the second collision of particle i with dF along (t(z, z', T), 
(r, 0, s, sin((/3)) - the coordinates of the collision, and 77^ is the projection to the 
coordinates of the z*'* particle. 
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We will prove Lemma [S] as well as all other technical lemmas in this argument 
at the end of the subsection. 

In the following, we will slightly abuse the notation by using the operator V* 
to push forward measures projected on components of D associated with a single 
particle or a particle and the disk. We operation is well defined since we will only 
do so on time intervals for which the chosen component does not interact with 
the rest of the system. 

In Lemma E] we show that if we start with a uniform measure v with density 
1 in a fixed-size neighborhood of a collision point of ct(z,z',T), then it's push 
forward at the next collision point has a uniform lower bound on the density in 
some uniformly sized neighborhood. 

Lemma 6. Let I be a line segment connecting two points r and r' in dF forming 
an angle ip with dF, satisfying a + e < \ sin((/j)| < 1 — e. Let s, Smin ^ s < Smax, 
be the speed drawn at r and let t he the time to trace the segment I with speed s. 
Let s' , Sniin < s < Smax, and if' , a + e < \ sin((/3)| < 1 — e, be the new speed and 
angle drawn at r' . 

Then there exist 77', C' > such that Vl^iH^(rfi,s.sin(-4>)) > mMH^j (r',o,s',sin(<p')) • 

It remains to treat the special case of a particle that sets the disk's angular 
velocity. 

Let To be the point of the second collision with dF of the particle chosen 
to set the disk. Let sq and sin((po) be parameters required for the particle to 
reset the disk to uj' . Let ri G dD be the location of the disk collision, ti be 
the time of flight between tq and ri , r2 € dF be the location of the following 
collision with dF, t2 the time of flight from ri to r2, and 9' be the position of 
the marked point on the disk at time ti + 12- Then the outgoing speed at ri is 
si = + Sq cos^ {ip'f^ ) — y^oj'^ + Sq — Sq sin((y9o)/ a and the outgoing angle ipi 
satisfies sin(95i) = lu/si. Denote the speed and the angle drawn at r2 by S2 and 

Lemma 7. Assume |a; — uj'\ > n for some n > and that sin(93o) < . Let 

Fluj' = — Ci ^' + Cui) X ~ Cui,^' + Clj) o,nd let hh^, &e the uniform measure 
on H^i with density 1. Then there exists r],^ > such that, if a particle resets lo 

to uj' , then 77a;[7'^+*^(MHa'~o:0,so,sin(i^o)) ^ he,'^))] ^ VuiVofJ-H^- Here LI^ denotes 
the projection to the disk coordinates. 

When a particle collides with the disk that rotates at a set angular velocity w, 
it 'looses' its density to the disk, i.e., LI,[Vl''^*^{^iHc^(rofi,so.sin(vo)) ^ ^{e^Lu))] is not 
absolutely continuous. Therefore, we need to ensure that the particle 'restores' 
its density once it collides with dF again. 

Lemma 8. There exists k > 0, 772 > 0, and C2 > such that if k ^ [^2k, 2k] and 
the particle is sent to change the angular velocity of the disk from uj to k or from 
k to Lo' , then, upon return to dF at time ti + 12, F[i[Vl^~^*'^{^H({ro.o,so.sin(<po)) ^ 

S{e,uj))] > »?2»?0Mfl'c2c('-2,0,S2,sin(ip2))- 

Combining Lemmas [5][H] completes the proof of Prop.[T] □ 
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Proof of Lemma [5} acquiring density 

Originally, each particle is assigned a point measure that evolves deterministi- 
cally until collision with OF at point r and then gets perturbed in s and sin((/3) di- 
rections. Then the measure supported on a two-dimensional sub- manifold evolves 
with billiard dynamics until next collision, when it gets perturbed in s' and 
sin((p') directions. Overall, the final coordinates (r', s' sin(<p')) depend on the 
original r as well as randomly drawn s, sin((p), s', and sin((^'). The Jacobian of 
this mapping does not depend on r and is equal to 

J _ dr' _ Ijv) Ijip) _ l{vf 

dsm{ip) ds cos2((y9) s cos^ [1^)3 

where l{if) is the distance of flight between collisions. This determinant is clearly 
bounded below and above by some J„iin and Jmax if Smin — C 1^ s < Smax + C 
and a + t — C < \ sin((/j)| < 1 — e + C,. Therefore the density in a neighborhood 
H(^{r' , 0, s', sin(iy9')) is bounded below by 770 := p^in/'^max for small enough C,. □ 

Proof of Lemma [6} Pushing density forward 

Let us first push the measure /^//^(r,o,s,sin(i^)) forward for time t under the 
billiard flow J-^ on the circle, i.e., the deterministic dynamics with the an- 
gles of incidence equal to the angles of reflection. Since the speed and the an- 
gle are preserved under the billiard flow, the Jacobian of in the variables 
(r, ^, s, sin(iy9)) is equal to 1. Therefore, the density of -7^i/^ff^(r,o,s,sin(v)) equal 
to 1 in the neighborhood J^*-ff^(r, 0, s, sin((/3)), which is 'skewed' in r and ^ vari- 
ables. Due to bounds on s and sin(9j), however, there exists Ci > such that 
i?Ci.c(r',0,s,sin((p)) C J"*(7I^(r, 0, s, sin((^)). 

Pushing the measure A'fl',(r,o,s,sin(i/3)) forward under 7^* keeping track only of 
the part that stays in a small neighborhood of (r', 0, s', sin((p')) is equivalent 
to pushing ^JiH■r,{rfl.s,sln{^p)) forward first by (with additional reflections off 
Lu U Ld, which only changes sign of the Jacobian), then perturbing in s and (f 
variables, and then applying the change of the chart map from a neighborhood 
of (r', 0, s, sin((p)) to a neighborhood of (r', 0, s', sm{if')). Indeed, the variables s 
and sin((^) are preserved under the billiard flow and we can think that the 
perturbation in s and tp variables occurs precisely at time t (and not necessar- 
ily at collision). Then we apply the change of the chart map T from a neigh- 
borhood of (r'j 0, s, sin((p)) to a neighborhood of {r' ,0, s' ,sm{(p')), which maps 
(r, ^, s, sin((^)) — >• (r', s', sin(</?')) as follows: f' = f — (r' — r), s' ~ s — {s' — s), 
sm((f') = sm{(f) — (sin((^') — sin((^)), and ^' = (1 + ■^^-y^)Ci where the last 
relation follows from It = S. 

s' s 

When we perturb ^//^^ ,.(r-',o,s,sin(v3)) in variables s and ip, the fraction of 
the measure that stays in H(^^.(^{r' ,0, s,sm{(f)) is 2pi„in(CiC)^, where pmin = 
min{ ^^ Sj^ij^e^'^'^mi", ^max^"^''"""'}- The Jacobian for the change of the 

chart mapping T is (1 -f- -^^-7-^), which achieves its minimum at |iais-. In addi- 
tion, Hf^^^^{r',0,s',sm{ip')) C THt^,.^{r' ,0, s,sm{ip)). Therefore,"'"' 

/ M//^{r,0,s,sin{(/;')) — ^Pmin^US. j ^Hs^^ Cr',0,s',sin((/?')) 

Lemma [HI follows. □ 
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Proof of Lemma [7) acquiring density for the disk. After collision, the 
disk variables satisfy lo' = sq sin((/3o)/a and 6' — 6q + wicoi + — tco\) = 
9o + oj't + {(jj — a;')icoi, where icoi = is the collision time of the particle and 
t = ti + 12. An easy computation shows that the Jacobian for this mapping is 

Juj = l[ipa)[cos[ipQ) - sm(<y9o)^^ 

aso otpo 



as 



-l{ipo)[cos{ipo) - sin(<y9o)tan((^o)]. 



If we choose | sin((/3o)| < '^'{^^1 ^ then there are uniform lower {J^i^) and upper 
(•^max) bounds on the Jacobian and there exists Cuj such that the pushed forward 
measure has density at least l/o/^,'^ax i'^ ^ neighborhood {9' — Cuj,0' + (^) x {oj' — 

Proof of Lemma [8} re-acquiring density for the particle that hit 
the disk. Similarly the Jacobian for the mapping from (s, sin((^), s', siii{(p')) to 
{r s' ,sm{(p')) can be computed given the particle resets the disk from uj to 
uj' . If we plug either w = or w' = 0, the Jacobian is bounded above and below 
by some positive constants. Therefore, there exists k > such the Jacobian for 
resetting to/from any k e [—2k, 2k] will also have lower and upper bounds. By 
the same reasoning as in Lemma[5l we conclude that there exist 772 > 0, and C2 > 

such that -n"i['P*i+*= (^H^(ro,0,so,sin(vo)) ^ ^{0,ij))] > V2ll0f^H^^^{r2fi,S2,sm(v2))- ^ 



3.5. Mixing. Mixing and convergence of initial distributions to the invariant 
measure follow almost immediately from our existence argument for the invariant 
probability measure. In the proof of Prop. [1] a lot of effort has been devoted to 
guarantee lower bounds on the densities of the pushed forward measures. A 
weaker property of Markov processes, irreducibility, can be shown in a similar 
manner by dropping the lower bounds on times and densities and allowing the 
paths we follow to start anywhere in the phase space. 

Definition 2. A continuous-time Markov process <l>t is called irreducible if for 
all z G i7, whenever Leb{A) > 0, there exists some t > 0, possibly dependent on 
both z and A, such that V*{z,A) > 0. 

Lemma 9. Markov process <l>t is irreducible. 

The proof of Lemma [H] is a simple modification of the proof of Prop.[T] 

Definition 3. The Markov process (!>t is called ergodic if an invariant probability 
measure v exists and 

lim \\V\z,,-) - v\\ =0, Vz e /2 

The proof of ergodicity of relies on sampling the Markov process <l't at 
integer times, which generates a discrete-time skeleton chain Denote the 
transition probability kernel for ^1 by Vi. The following theorem by Meyn and 
Tweedie relates skeleton chains to the ergodicity of the Markov processes. 
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Theorem 3. f21[ Thm. 6.1] Suppose the Markov process <Pt is irreducible and v 
is an invariant probability measure for <Pf Then (Pt is ergodic if and only if <Pi 
is irreducible. 

Proposition 2. The Markov process (!>t is ergodic. 

The mixing of the invariant measure for the Markov process and the con- 
vergence of initial distributions to the invariant measure follow from ergodicity 
by the Dominated Convergence Theorem. Indeed, 



and to show mixing one may replace Q with B in the second integral. 
Proof of Prop. [H 

Ergodicity follows once we show that the time-1 sampled chain is irre- 
ducible. The proof is a modification of the proof of Prop. [TJ What we really need 
to show is that for any z, z' G i7, there exists T G N and a sample path (t(z, z', T) 
from z to z' in time T . In addition, some density is acquired and carried through 
along the path. No lower bounds are required. The existence of a sample path 
is guaranteed in the proof of Lemma |4] noting that allowing z and z' to be not 
in C only changes times it takes for particles to reach dF for the first time in 
forward or backward times respectively. Time T for Lemma [3] is allowed to vary 
in a range of values. In particular, we can choose T such that T is integer (here 
is does not have to be the same for all states either) . Modifying Lemmas [5]|8] to 
apply for all injection parameters and dropping the lower bounds, we conclude 
that (Pi is irreducible. □ 

4. Sub-exponential mixing 

Proposition 3. There exist (many) initial probability distributions X on Q that 
converge to the unique invariant measure fj, with sub- exponential rates, i.e. 3"; > 
such that for T large enough 



In particular, the unique invariant measure /i for the Markov Process <Pt is not 
exponentially mixing, i.e. there exist (many) Borel sets A, B such that 



Proof of Prop. [3l 

Let Bt be a set of z g J7 such that at least one particle in z does not have 
any collisions for time T. We are interested in estimating ^{Bt). 

The proof is rather similar to the sub-exponential mixing proof for a system 
driven by thermostats in [23 except the dynamics on Bt is not deterministic 





B 
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and there is no potential to aid the estimates on an upper bound on the measure 
ofBr- 

Prop. [T] ensures that there exist To,r] > such that 'P'^°{x,-) > rpnc, where 
> and mc is the uniform probability measure on C. By the invariance of /i, 
we conclude /i > ^{C)rjmc ■ 

Let 5 < d/ (smax ■ This guarantees that for any z e C, any particle 

in z that experiences a collision with dF at time t < S does not experience any 
disk collisions on time interval [0,t). 

Denote by f2i the projection of the phase space f2 into components associated 
with i*'' particle and let ttiq. be the uniform probability measure on Sli. Let 
Gi = {{ri,Si,ifi,^i) : the particle will collide with dF in time S}. Here we use 
coordinates defined in subsect. 13.21 Let 7 = Qi{Gi) > 0. Then if z is drawn 
uniformly in C, with probability 1 — (1 — 7)*^, at least one particle in z will hit 
dF in time S. 

We are interested in the probability that a randomly emitted particles will 
not collide with the disk or dF in time F + 6 > F + {S — t). Consider only the 
situations when a particle will fly at least distance d before the next collision, i.e. 

I sm{(p) I < y^l — (-^^)^, which happens with probability y^l — i^^Y- Then 
the probability that a randomly emitted particles will not collide with the disk 
or dF in time F + 5 is greater or equal than 




R + d Jq 

T--fErf[^].l^^^ + 0(J-)>^ 



Then at least (1 — (1 — 7)*^) x ^-fraction of mc, and therefore at least 



for some a > and F large enough. 

lx{C)rj{l— (1 — 7)'^) X ^ fraction of ^, ends up in Bt in time 5. By the invariance 
of /i, we conclude that ij.(Bt) > /i(C)7](l — (1 — 7)*^) x ^ for large enough F. 



To estimate an upper bound on the fraction of /x that will end up in Bt, 
we observe that the only way to get to Bt not from Bt is to emit a slow 
enough particle from dF. Therefore, starting from any initial distribution A, the 
probability to end up in Bt is bounded above by 



4/33/2 



s'e'^^'ds < ^ 



for some a' > and F large enough. 

Thus we obtain that for any n G N and F large enough, h{Bt) — ^J,{BT+ns) + 
nf{F), where /(T) « ^. Here / « 5 means /(T) = 9{9{F), i.e. there exists 
> such that e'ff(r) < f{F) < ^g{F). 

00 

Summing up, we conclude that ^{Bt) ~ {T+kS)^ ~ W- 

k—O 

Note that the dynamics for our system is statistically very similar to the 
dynamics of an expanding map with a neutral fixed point (aka the Pomeau- 
Manneville map). Indeed, the mass originally in Bt evolves in Bt for at least 
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time T; extra mass is 'deposited' from the parts of the phase space where particles 
experience nearly tangential collisions. One can complete the proof of Prop. [3] 
using an argument very similar to [23 and [24 . 

Let X ^ fi with dX — ipdfi be such that (p > 1 -f- c on Bto for some small 
c> 0. Then for kS > Tq 



IIKA - Mil > (K^KBks) - ^KBk5) > XiBks+ns) - KBks) 

> (1 + c)^iiBks+ns) - KBks) = [(1 + c)^^i^M±!i^ _ iHBks) 

r/1 I ^ KBkS+ns) 11 /D \ 

f4Bks+nS) + nfn[ko) 

Then for N large enough and k = Nn 

IIKA - Mil > [c - (1 + c)—-^^^^^^^l——MB^,,s) > " ^' 



□ 
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